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1. INTRODUCTION 



Over the past few years, the interest in littoral, shallow water environments has 
increased significantly The success of many naval systems in these environments depends 
on our ability to predict and interpret shallow water acoustic propagation The sound 
speed profile in shallow water typically consists of an upper mixed layer below which the 
sound speed decreases to the bottom Typical shallow-water environments are found on 
the continental shelf with approximate depths of 200 m. The wave propagation is 
downward refracting, that is, bottom-interacting paths dominate in long-range propagation 
as depicted in Figures 1 and Figure 2. Etter (1991) summarized the effect of the sea floor 
as: 

• forward scattering and reflection loss (but complicated by refraction in the 
bottom), 

• interference and frequency effects; 

• attenuation by sediments; 

• noise generation at lower frequencies due to seismic activity; 

• backscattering and bottom reverberation 

The purpose of this thesis is to quantify the influence of bottom variability on the 
resolution of pulse transmissions in simple, shallow water environments. 

Rough bottom interfaces can be modeled by assuming some realistic form of the 
spectral shape. The forward scatter of the acoustic field can then be treated 
deterministically Fluctuations of the sediment properties, such as porosity, density, 
permeability, and shear strength, determine the sound scattering from the bottom volume 
Yamamoto (1995) summarized the relationship between spectral strength and sediment 
type. We used this information to generate random two-dimensional wave-number 
spectra of the velocity variability in order to simulate a realistic ocean sediment volume. 
These techniques are described in more detail in Chapter II 



Typical sound speed profile in shallow water 




Sound speed (m/s) 



Figure 1. Typical sound speed profile in shallow water. 



Ray trace in shallow water 




Figure 2. Ray trace in shallow water. 
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Wave propagation in an inhomogeneous medium is defined by the wave equation 
The solution can be obtained by several approximate methods, for example, Wentzel- 
Kramers-Brillouin (WKB) approximation, ray approximation or parabolic equation 
approximation. The last method was used throughout this thesis study. Over the past few 
decades, ocean acoustic models have been developed and improved to analyze measured 
data and simulate sound propagation in real ocean environments Modeling provides an 
efficient means by which to parametrically investigate the projector (source) radiated 
sound pattern under the influence of environmental conditions A perfect model is one that 
perfectly represents reality In practice, such a perfect model would be as complex as the 
problem it is attempting to represent So, there are some assumptions built into every 
computer model This thesis uses a version of the PE model developed at the University of 
Miami known as UMPE (Smith and Tappert, 1994) 

In order to investigate the effects of bottom interface roughness and volume 
fluctuations, we simulate broadband pulse propagation for four different scenarios as 
follows. 

• flat bottom with homogeneous sediment; 

• rough bottom interface with homogeneous sediment; 

• flat bottom interface with realistic sediment sound speed fluctuations, and 

• rough bottom interface with realistic sediment sound speed fluctuations 






4 



n. BACKGROUND 



A. THE ACOUSTIC MODEL 

Because propagation ranges are typically much larger than ocean depths, it is 
appropriate to use a cylindrical coordinate system for long range propagation Customarily 
we define j as depth, r as range outward from the origin, and tf> as the azimuthal angle 
{()—>2n) defining direction along a radial in the horizontal plane as depicted in Figure 3 In 
general, full 3-D modeling of the propagation is seldom done due to mathematical and 
computational complexity. It is often a good approximation to the real ocean acoustic 
environment to consider it azimuthallv invariant, i e , to assume that environmental 
parameters such as sound-speed profile, water depth, and bottom composition are 
independent of azimuthal angle to the source. With this assumption, we concentrate the 
analysis in depth versus range. 

As mentioned in the introduction, the PE method is used in this thesis to model 
acoustic propagation. The parabolic equation method was introduced into underwater 
acoustics in the early 1970s by Hardin and Tappert (1973) who devised an efficient 
numerical solution scheme based on fast Fourier transforms We will derive a PE solution 
in a 2-D constant-density fluid environment where we have assumed azimuthal symmetry 
and hence no dependence on the <f> coordinate. 

We begin with the wave equation for a constant-density medium in cylindrical 
coordinates (r,z) independent of <j> 




0 ) 



We now assume we have a steady state, time-harmonic acoustic field 

p(r,z,mt) = p(r,z)e~' m . 



( 2 ) 
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Figure 3. Cylindrical coordinates in ocean. 



6 



( 3 ) 



Substituting this into Eq. (1) leads to the Helmholtz equation, 

d 1 p 1 dp d 2 p 2 i 

— f + + — f + ^ 0 V / 7 = ° , 

dr~ r dr dz 



co 



where p(r,z) is the acoustic pressure. Ay = — is the reference wave number, 



//(/', r) = t is the index of refraction, c 0 is the reference sound speed and c(r,z) is the 

Ma-J 

acoustic sound speed It is within c(r,z) that all features of the environment are 
represented 

Tappert (1977) assumed the solution of Eq (3) in terms of an outgoing Hankel 
function with an envelope function y/(r,z) 

p{/\z) = y/{i\z)H\;\k {) r) , (4) 



where y/(r,z) is assumed to be slowly varying in range. In the far-field (k Q r» 1), the 
Hankel function can be approximated by 







( 71 \ 

V» r -7J 



(5) 



thus an alternative relationship between the acoustic pressure, p(r,z), and the slowly 
modulating envelope function, y/(r,z ), is 

P(dz) = P^-yy/ir^zje'^ . ( 6 ) 



This is the standard definition of the so-called "PE field function" y/ scaled such that at 
/' = R , t , \y/\ - 1 and |/r| = P n . Substituting (6) into (3) yields 



o_ 

r y/ 
cr 






(7) 



Now we introduce the crucial paraxial approximation based on the slowly modulating 
character of y/. 



d 2 y/ 

~dd 



« 2 ik, 



dy/ 

dr 



( 8 ) 
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We then obtain the following wave equation. 



2/A 



ri// c~y/ 



— -r A,:(/r - l)y/ - 0 



cr oz 

which is the standard parabolic equation This can be rewritten as 



: , _Lif_ 

» 2 k" k; A 1 






( 9 ) 



( 10 ) 



which is the form first introduced to the underwater acoustics community by Tapper! 
(1974) 

Equation (10) is a one-way wave equation which can be solved by a range- 
marching numerical technique. The UMPE model applied in this thesis uses the split-step 
Fourier (SSF) method (Hardin and Tappert, 1973) which solves the parabolic equation by 
Fourier transform techniques. Define the complex Fourier transform pair as 



(//(/-. r) = j k : )e‘ k: 'dk : = h '[y/(/\A\)j , 

where k : is the vertical wave number We may rewrite Eq. (10) as 




(H) 

( 12 ) 



( 13 ) 



where 



7op 2k: d z 2 



( 14 ) 



represents the kinetic energy operator and 

■ 05 ) 

represents the potential energy operator This representation of the operators as kinetic 
and potential energy operators is especially insightful when one wishes to form the ray 
equations which have I lamiltonian form (Smith, 1992) 

In the PE/SSF method, solutions of y/are marched out in range by 
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y/(r + A r) = [f/ ( r ) (16) 

The operator (J op is simply a multiplication operator in z-space but T op contains cross 
terms which couple different depths However, T „ p is a multiplication operator in vertical 
wave number space It is therefore desirable to separate the application of each operator, 
one in r-space and one in A'-space (here, k is the vertical wave number) In other words, 
the PE field function y/ is specified at some range r initially in the r-domain A 
transformation is made to the A'-domain followed by a multiplication of the A-space 

operator c ,< "‘ v 7 "'' where T 0f , is the corresponding kinetic energy operator in the wave 
number domain The result is then transformed again to the .--domain and is followed hv a 
multiplication of the r-space operator c ' k '"' The final result is the field function at 
/• Ai\ that is, 

y/(r + A r,z) = x F~ ] ^^ rJc) x [F(y/(/-,z))][ . (17) 



This is the general algorithm of the PE/SSF method. 

Previously, we have assumed the operators took the forms defined by Eq (14) and 
( 15) These forms are commonly referred to as the “standard PE” and are only one set of a 
number of various operator forms In our simulation, we will use a higher order form 
which was introduced by Thomson and Chapman (1983) referred to as the “wide-angle” 
PE approximation (WAPE) with operators given by 



L 




( 18 ) 






(19) 



We then need to define the initial conditions for the PE field function, y/fr 0,z) It 
is desirable to begin the calculation by specifying the source in the Ar-domain Thomson 
and Bohun (1988) have shown that a wide angle source can be defined as 
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( 20 ) 



1 f o p ^ 

i//(0,k) = Ak — I 






I I 



2 V Jtk n ) \ v 

where A- k 0 corresponds to a propagation angle of 90° 

B. TRAVEL TIME ANALYSIS 



UJ 



(l*M„) - 



I he calculation of the time domain arrival structure at a given range is an integral 
part of the analysis in this thesis The pressure field in the time domain is represented by 

p{r,zj) = F-\p m {r,z)\ = Y. a PJ< r ’ z y a ’ t > ( 21 ) 

where p a (r,z) is a single frequency component. The arrival time structure at some fixed 
range r=R can be computed using the complex field p a (R,z) for many frequencies and 

then Fourier transforming to obtain p{ R,z,t), the set of complex pressure values in the 
time domain To avoid large transform sizes, the UMPE model basebands the solutions 
and shifts the center frequency bin to zero To reduce the effect of side-lobes, we apply a 
Hanning filter over the bandwidth before transforming to the time domain Additionally, 

K 

ICO — 

we neglect the phase factor e' k ‘ ,H - e e * by simply performing the Fourier analysis on the 
field function i//„/R,:). Therefore, the time domain is heterodyned around the value 
l„ R c„ and arrival times are computed as values of “reduced time” or (/-/„). 

C. ROUGH BOTTOM INTERFACE 

The ocean bottom roughness is assumed to be characterized by a 2-D power law 
spectrum at high wave numbers which can be written as 

( i ^ 

W 2 \k »- — J = ak~ p , (22) 

where L a „, is the correlation length of the roughness and /? is the spectral exponent (Note 
that our use of wave number is now with respect to range scales of interface roughness ) 
To give the full spectrum a realistically smooth structure, we assume the spectral form 
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( 23 ) 



K(k) = ~ T 

(l + 2 ) : 

The nns roughness represents the average amplitude of the undulated bottom interface 
variation and the correlation length can be thought of as a scale of self-similarity. The 
normalization factor /j is defined in terms of the rms roughness a, i.e., 

<jr • 

27tJ W 2 {k)kdk - a 2 (24) 

(i 

and therefore 

I />' 

// -• - - - I ct / ( 25 ) 

n 2 

We only need the 1-D roughness spectrum along the slice in the two dimensional 
profiles (depth z and range r). Along the direction of range propagation in cylindrical 
coordinates, we have 

wtf) = -j J==_ =,_ W,(K)dK , ( 26 ) 

where K is the total horizontal wave number and -A' is the wave number along the 
propagation direction. Substituting Eq (23) and Eq (25) into Eq (26), we obtain 

»;(*)= 4A) = • ' • (27) 

where 




and r(x) is the gamma function 

To avoid the annoying gamma function, the UMPE model assumes the 1-D 
spectrum has the form 

»;’(*)=( • 



1 1 



(29) 



The roughness realization is then simply rescaled by the desired value of a to provide the 
appropriate rms value 



I). SEDIMENT COMPOSITION 

Yamamoto (1995) summarizes the transformation between three-dimensional 
wave number spectra and the velocity variabilities which determine the ocean sediment 
properties This transformation can be simplified for our 2-D case. The 2-D congressional 
wave velocity staicture in the sediment may be written as 



where c(x,z) is the velocity at a point (x,z), c 0 is the average velocity at the seafloor and 
Scfx.z) is the relative variability. The relative variability is modeled as a zero-mean 
quasistationary random function, i.e.. 



The 2-D power spectrum of the relative velocity variability is denoted by S(k^k y j By 



It is a plausible assumption that the vertical power spectrum may be given by a power-law 
form 




(30) 




definition, the variance (S ]) represents the total energy of the power spectrum 




(33) 




where /? is the spectral exponent and B is the spectral strength constant. The 2-D power 
spectrum given by Yamamoto is 




(35) 



where a is the aspect ratio and /? 2 must be an integer 
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To produce a random realization of sediment volume fluctuations, we assume the 
complex wave number spectaim has a random amplitude and phase. Ac' 9 ' for each wave 
number component. The random amplitude of each component of the power spectaim can 
then be obtained by A" =-ln(r) where r is an independent, uniformly distributed random 
variable in the interval (0,1). The random phase (f> is another independent, uniformly 
distributed random variable in the interval (0,2k). To avoid infinite magnitudes at low 
wave numbers, the spectrum is high-pass filtered The filtered spectrum is then Fourier 
transformed to obtain a realization of the sound speed variability The entire field of 
sediment fluctuations is first computed on a fine grid and then interpolated to match tite 
computational grid of the PE/SSF algorithm The entire scheme can be described as 
follows 



Step 1 . Construct the 2-D wave number field (k r (//), k, (/?/)) , according to 



M = 



in 



(36) 



(.(>/) = (//- i)aa. 



w = 1,2 vV . 



A, 271 

M< : = 7 7 , 

l 2 ma (37) 

k. ( m ) = (m - 1) A k z m = 1,2, . . . , M - M bol , 

where 

/• max = maximum propagation range, 

N = FFT size in range direction, 
r - maximum compution depth. 

: luil average bottom depth 
M = FFT size in depth direction, 

M hiil = index of average bottom interface 



Step 2. Generate two random numbers series for the random phase and 
amplitude, 

A{n,m) = -J- ln(a n m ) , (38) 
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(39) 



4(n,m) = 2nh nm 



where 

a,b - unifonn random number distribution 0 — > 1 






Step 3. Compute the random phase and amplitude for every mesh point in the 
2-D realized wave number field, S'(k (n)k .(///)) , i e , 

.(//), k (///)) - S{l<, (ii),k_(iiif)A{u,ni)e'* { ' u "' . (40) 



• Step 4. Add high-pass filter to normalize the exponent amplitude at low wave 
number 

• Step 5. Compute 2-D FFT to yield relative velocity variability, 

S c (r,z) = 2DFF1^S'(k r (n)k z {m ))] . 

The velocity structure of the sediment realization is then given by: 

c(/\ z) = c„ + c n X R.e(<5(/\ z)) . (41) 



E. BEAMFORMING 

In underwater acoustic propagation, much of the information about the 
environment and the relative source/receiver locations is contained in the vertical structure 
of the acoustic field. This is especially true in shallow water environments where the 
vertical structure can be defined in terms of a relatively small number of modes of the 
acoustic waveguide Because these modes are generally dispersive, the ability to resolve 
the vertical structure of the time arrival at the receiver can lead to estimations of the 
corresponding source location. 

Compared to the directivity of single sensors, a beamformer permits one to listen 
preferentially to wavefronts propagating from one direction over another. In this sense, a 
beamformer implements a spatial filter (see, for example, Defotta el al ., 1988, Ziomek, 
1995) Consider a line array of equally spaced elements positioned along the y-axis in an 
infinite, homogeneous medium with an undetermined number of remote sources Let the 
output of an element located at the origin of coordinates due to the llh source be s/t). 



14 



Under the assumption of plane-wave propagation in the far-field, a source from direction 
0 , produces the following sensor outputs 



'»(/) - s i 



//c/sin^)^ 



( 42 ) 



where // is the array element index, d is the spacing, G t is the angle of arrival from the Ith 

source, and c is the speed of propagation Beamforming is accomplished by applying 
weights and time delays to the individual element signals, and then coherently summing 
them The beamformer output when viewed in the direction is then 

c/(sin(6 , ; ) - si n(<9„, )) 



h ,S>) = Z M V V / 



nu 



( 43 ) 



where the ir n are weights applied to each element’s signal 

The directional response characteristics of the beamformer can be derived by 
considering a complex, plane-wave cw signal propagating across the array, 

•V(') = (44) 

The source produces the following sensor outputs due to the arrival angle 



>ico f t inkidsm(0 e ) 



«.(0 = 



( 45 ) 



2 n co e , 

where k f = — - = — is the wave number of the 1th source. The beamformer output is then 

Af C 



bJJ) - w n e mk, ^' n{l,r) s ' n(Wm, l - 



( 46 ) 



where 



.v 1 

;nk ,d Isinf 0, i sinf )j 

u e 1 

n 

/?— n 




( 47 ) 



represents the spatial Fourier transform of the array weights 

A two-dimensional FFT (spatial and temporal) can be used to implement the 
beamforming operation for an equally spaced line array. For example, consider a line array 
having N 8 equally spaced elements as depicted in Figure 4 The signals from each 
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element are sampled simultaneously at rate f s to obtain a multichannel time series x nm . M- 
point FFTs are taken along each channel to obtain the temporal frequency response, X„(j), 
where / is the frequency bin number. 1 e , 



•V„( /) = T „c , 0 < / < M - I (48) 

0 

The A '„(/) are represented in the figure as N horizontal rows of M FFT bins. 

Next, the N cross-channel FFT is performed with appropriate weighting along the 
array, i.e., 

x{j,k) = y w n XXjY‘ [2%]nk ,0<k<N-\ (49) 

n=l) 

The result of the N cross-channel FFT series to spatial frequency (or wave number) 
provides the phase difference information in depth for specific temporal frequencies. The 
relationships between the various indices and the variables used in the beamformer output 
are 



./, 



XT' 

k 2 n 
N d 



f 



6>.„ = sin"' 



c k , 



f { ) 



(50) 

(?») 

(52) 



where j and k are the indices in the frequency and wave number domains, respectively. In 
practice, 6 m is determined by interpolating between values of ki . After the data has been 
interpolated to arrival angle/frequency space, an inverse FFT is taken along each angle bin 
to provide arrival angle as a function of travel time 
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Figure 4. Beamforming example. 
From [DeFotta, 1988] 
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III. MODEL SIMULATIONS 



As was described in the introduction, the principal characteristic of most shallow - 
watei propagation is dominated by bottom-interacting paths Sound transmission in 
shallow water has an optimum frequency of propagation at longer ranges as shown by 
Jensen and Kuperman (1983) Optimum frequency is a general feature of ducted 
propagation in the ocean It occurs as a result of competing propagation and attenuation 
mechanisms at high and low frequencies. In the high-frequency regime there is increasing 
volume and scattering loss with increasing frequency. At lower frequencies the situation is 
more complicated. With increasing wavelength the efficiency of the waveguide to confine 
sound decreases (the cutoff phenomenon). Hence propagation and attenuation 
mechanisms outside the waveguide (in the seabed) become important. In fact, the 
increased penetration of sound into a lossy seabed with decreasing frequency causes the 
overall attenuation of waterborne sound to increase with decreasing frequency. Thus we 
get high attenuation at both high and low frequencies, while intermediate frequencies have 
the lowest attenuation As Jensen's study shows, the optimum frequency is in the range 
200-800 Hz for a water depth of 200 in and we will use 200 Hz as our center frequency 

We assume that the ocean bottom sediment is a fluid-like half-space which means 
that it supports only one type of sound wave, i e , compressional waves The information 
required for model simulations should include the following depth-dependent material 
properties: the sound speed profile in the water column and compressional wave speed, 
attenuation and the density in the seabed. For all cases, we will assume the water column 
is homogeneous with a constant sound speed. 

This chapter presents the model simulation results of four test cases. Case 1 
examines the propagation of a broadband, impulsive source over a flat, homogeneous 
bottom Case 2 examines the propagation in a similar environment with a rough bottom 
interface We then examine which parameters, i.e., rms roughness, correlation length or 
ratio of them, dominates the beam resolution Case 3 is also similar to Case I but with 
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realistic bottom volume lluctuations The factois which most strongly mtluence resolution 
m Case 2 and C'ase 3 will be combined into Case 4 

The basic range-independent parameters used in the UMPE model are listed in 

Table 1 



Model Type 


Wide Angle PE 


Source Function Type 


Wide Angle Source 


Reference Sound Speed 


1500.00 m/s 


Bottom Depth 


200.0 m 


Maximum Computational Depth 


400.0 m 


Mesh Points For Positive Depths 


512 


Maximum Range 


10 o km 


Step In Range 


2.0 m 


Mesli Points In Range 


5000 


Source Depth 


100.0 m 


Center Frequency 


200.0 Hz 


Frequency Bandwidth 


100.0 Hz 


Number Of Frequencies 


256 


Sound Speed In Water 


1500 m/s 


Sound Speed In Sediment 


1600 m/s 


Densit\ In Sediment 


1 500 kg/m ’ 


Attenuation In Sediment 


0 2 dB/m/kHz 



Table 1. IIMPE Model Parameters 



A. CASE 1: RANGE INDEPENDENT 

We begin with the simple case consisting of a constant speed water column and 
fluid-like homogeneous sediment separated by a flat interface The corresponding 
transmission loss field for the center frequency is depicted in Figure 5 We note the 
symmetric constructive interference due to the finite number of trapped, propagating 

normal modes By Snell's law, there is a critical angle at arccos— — =20 4° (relative to 

1600 

horizontal) below which no energy is transmitted into the bottom. From the UMPE output 
field with beamforming at 10,000 m, we obtain the angle-time plot shown in Figure 6. As 
expected, this energy arrives primarily near horizontal with an aperture of approximately 
±20° The steeper-angle energy >20° is rapidly lost by transmission into the seabed. 
Furthermore, the lowest modes arrive first, before the higher modes, as would be expected 
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Case 1 (Flat bottom) 
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Figure 5. Range-depth transmission loss field for flat bottom interface. 
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Figure 6. Arrival angle-time plot for flat bottom interface. 
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from an analysis of the modal dispersion in this environment A standard calculation of the 
number of trapped modes was found to be consistent with the number of arriving modal 
peaks when the range of the transmission loss palette was extended to expose the weakest 
modal arrivals. 

B. CASE 2: ROUGH BOTTOM INTERFACE 

A rough seafloor will cause attenuation of the mean acoustic field propagating in 
the ocean waveguide due to incoherent scattering thereby decreasing the mean signal to 
noise ratio For individual pulse transmissions, the instantaneous signal becomes smeared 
in time and space due to mode coupling A rough bottom interface is added to the 
environment of Case 1 in order to investigate the effect of scattering on the instantaneous 
resolution of the pulse in the water column 

As described in Chapter II, there are two parameters that control the roughness, 
i e the rms roughness and the correlation length In order to determine the dominant 
factor we simulated two sets of comparisons In one, the correlation length is fixed at 
1000 m and the rms roughness has values of 2, 5, and 10 m. In the other, the rms 
roughness is fixed at 2 m with correlation lengths of 40, 200, and 400 m. Figure 7-12 
show the beamformed angle-time plots at 10,000 m. We can see that increasing the rms 
roughness and decreasing the correlation length will reduce the resolution but it is 
difficult to see which factor dominates However, the distortion changes look similar and 
suggest that the ratio of rms roughness and correlation length may be the important 
parameter We simulated several realizations with fixed ratio of the rms roughness and the 
correlation length to confirm our hypothesis, specifically, (rms roughness) (correlation) = 
I 500. I 100 and 1:20 The beamformed results are displayed in Figure 13-1 5 Combined 
with the previous simulations, we find that higher ratios retain higher resolution and lower 
ratios degrade the resolution more severely Figure 16-18 show single frequency (200 Hz) 
transmission loss fields in depth-range slices for various rough interface realizations 
Compared to the Hat bottom case, we see that there is more energy penetration into the 
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Figure 7. Arrival angle-time plot for rough bottom interface with rms roughness = 
2m and correlation length = 1,000m. 
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Figure 8. Arrival angle-time plot for rough bottom interface with rms roughness = 
5m and correlation length = 1,000m. 



24 



TL (dB) 



o 




o 

“ (N 

i 



o 


o 


o 


o 


o 


o 


o 


d 


o 


o 


o 


d 


o 


o 


05 


CD 


CO 




CO 


CD 


05 



(gap) ajguB ^BALuy 



Figure 9. Arrival angle-time plot for rough bottom interface with rms roughness = 
10m and correlation length = 1,000m. 
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Figure 10. Arrival angle-time plot for rough bottom interface with rms roughness = 
2ni and correlation length = 40m. 
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Figure 1 1. Arrival angle-time plot for rough bottom interface with rms roughness = 
2m and correlation length = 200m. 
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Figure 12. Arrival angle-time plot for rough bottom interface with rms roughness = 
2m and correlation length = 400m. 
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Figure 13. Arrival angle-time plot for rough bottom interface with rms roughness = 
5m and correlation length = 2,500m. 
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Case2 (rms=2m, correlation length=200m) 
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Figure 14. Arrival angle-time plot for rough bottom interface with rms roughness = 
2m and correlation length = 200m. 
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Figure 15. Arrival angle-time plot for rough bottom interface with rms roughness = 
5m and correlation length = 100m. 
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Case 2 (rms=5m, corr.=2,500m) 
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Figure 16. Range-depth transmission loss field for rough bottom interface with rms 
roughness = 5m and correlation length = 2,500m. 
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Case 2 (rms=2m, corr.=200m) 
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Figure 17. Range-depth transmission loss field for rough bottom interface with rms 
roughness = 2m and correlation length = 200m. 
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Case 2 (rms=5m, corr.=100m) 
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Figure 18. Range-depth transmission loss field for rough bottom interface with rms 
roughness = 5m and correlation length = 100m. 
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bottom at longer ranges, and the stable interference structure has broken down into a 
complex pattern. 

C. CASE 3: SEDIMENT VOLUME FLUCTUATIONS 

The scattering due to the velocity variability within the sediment is considered to 
be the major scattering mechanism in many geological fields in shallow water based on 
extensive measurements of acoustic backscattering and sediment cores (Jackson and 
Briggs 1992, Yamamoto 1995). The geology of shallow-water seabeds is complicated 
because of the physical, chemical, and biological activities which took place though the 
long geological time spans. The measurements of the velocity and density variabilities in 
the ocean are beyond this thesis subject. For our simulations, we shall use the sediment 
model described in Chapter II with some typical values extracted from Yamamoto's 
work. 

In this case, we use a flat bottom interface and add the sediment volume 
perturbation. According to Yamamoto's model, the bottom sediment perturbation is 
controlled by three parameters: spectral exponent (3, aspect ratio a, and the spectral 
strength constant B. In his conclusions, the spectral exponent P and aspect ratio a are 
fairly constant regardless of the sediment type but the spectral strength constant B varies 
depending on the sediment type and location. We will fix a and P and vary the values of 
B by choosing three typical spectral strength constants, e.g., 2x1 O' 3 , 2x1 O' 4 and 2x1 O' 5 
(units of B are /w 1 '* 3 ). The sediment Field for strength constant B = 2xl0‘ 3 is depicted 
in Figure 19. The rms value of the sound speed fluctuations in the sediment for this value 
of B is 22.4 m/s. The beamformed results for the three values of B are shown in Figure 
20-22. The pulse shape does appear to degrade slightly as the spectral strength constant is 
increased Overall, however, the sediment volume perturbation doesn't seem to affect the 
pulse resolution significantly. 
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Figure 19. Sediment field for strength constant B=2xl0 ' m 1 
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Amplitude of velocity variability (m/s) 
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Figure 20. Arrival angle-time plot for sediment volume fluctuation with strength 
constant B=2xl0' 3 m l p . 
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Figure 21. Arrival angle-time plot for sediment volume fluctuation with strength 
constant B=2xl0 4 ni H \ 
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Figure 22. Arrival angle-time plot for sediment volume fluctuation with strength 
constant B=2xlO r, m 1 p . 
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D. 



CASE 4: KOIJGH BOTTOM INTERFACE WITH SEDIMENT VOLUME 
FLUCTUATIONS 



We have seen the effect of rough bottom interfaces and sediment volume 
fluctuations separately. We now need to combine these two factors From Case 2 and 
Case 3 results, it appears that the effect of a rough bottom interface dominates the 
pulseresolution The strongest sediment volume perturbation strength constant, B = 2x10 '’ 
m 1 ' 1 ', was chosen to produce the largest volume effect. The interface correlation length 
was fixed at 1000 m and the rms roughness was varied from 2 m to 10 m. Comparing 
these results to those cases with no volume fluctuations showed little difference indicating 
that these levels of interface roughness dominate the effects on pulse resolution 

E. PULSE RESOLUTION ANALYSIS 

We ha\e shown that the effect of the bottom roughness and sediment volume 
fluctuation degrades low frequency pulse resolution in shallow water propagation. In this 
section, we will quantify the effect of rough bottom and sediment volume fluctuation on 
pulse resolution in order to obtain a clearer relationship. The conventional angle resolution 
is defined as an angle width of Me peak amplitude. In order to simplify the Case 2 analysis, 
we fixed the correlation length at a value of 1000 m and varied the rms roughness. Figure 
23 shows several slices in reduced times of 1, 1.5, 2, 2.5 and 3 seconds from an arrival 
angle-time plot for rms roughness of 0 (the Case 1 results), 2, 5, 10 and 20 m. Each 
reduced time data point is actually an average of five points spread over a time bin width 
of 0.2 sec, i.e., the datum for the reduced time of 1 sec is an average of data at 0.90, 0.95, 
1.00, 1.05, and 1.10 sec. Polynomial fits to the data produce the smooth curves shown in 
the figures We can see a good angle resolution at lower rms roughness values. As the rms 
roughness is increased, the resolution is found to degrade, particularly for the later arrival 
times This is expected since later arrivals correspond to higher modes which interact more 
stiongly with the bottom As an estimation of this effect, we note that the angular pulse 
width appears to double for an rms value of ~10 m corresponding to an rms: correlation 
ratio of 1 100. 
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In Figure 24. slices are taken at arrival angles 0, 5, 10, and 20 deg In this case, 
three points over a width of 1 deg are averaged to produce each data point shown. There 
is an anomalous point at an rms of 2 m due to unresolved peaks of two modes. This point 
was not used in the curve fit. The time resolution degrades when the rms roughness 
increases for all arrival angles except 0 deg. This is due to less bottom interactions of the 
horizontal propagation. Furthermore, we notice that the time width does not degrade as 
rapidly as the angular width None of the rms roughness values appear to double the pulse 
width in time I'he peak minimum transmission loss also increases with rms roughness as 
shown in Figure 25 by roughly 5-10 dB at an rms value of 10 m 

To perform a similar analysis of the Case 3 environment, the spectral exponent and 
aspect ratio were fixed Obviously, the sediment volume fluctuations show more stable 
results with varying strength constant as depicted in Figure 26-28. The angle and time 
resolution change little with increasing strength constant. The peak (minimum) 
transmission loss also remains stable. 

In Figure 29, the peak minimum transmission loss computed with the combined 
effects of a rough bottom interface and sediment volume fluctuations shows very similar 
results as in Figure 25. It is reasonable to conclude that these levels of bottom roughness 
dominate pulse resolution effects in simple shallow water propagation. 
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figure 23. Angle width versus rms roughness for reduced time =1, 1.5, 2. 2.5 ;ind 3 
second. 
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Figure 24. Time width versus rms roughness for arrival angle = 0, 5, 10, 15 and 20 
degree. 
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Figure 25. Peak minimum transmission loss versus rms roughness for reduced time 
= I, 1.5, 2, 2.5 and 3 second. 
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Figure 26. Angle width versus strength constant for reduced time = 1, 1.5, 
3 second. 
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Figure 27. l ime width versus strength constant for arrival angle = 0, 5, 10, 15 and 
20 degree. 
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Figure 28. Peak minimum transmission loss versus strength constant for reduced 
time = 1, 1.5, 2, 2.5 and 3 second. 
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Figure 29. I’eak iiiiniinuiii transmission loss versus rms roughness for reduced time 
= I. 1.5, 2, 2.5 and 3 second combine with strength constant = 2x1 O' 3 m 1 . 
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IV. CONCLUSION 



In this thesis we have presented a numerical investigation of the effects of bottom 
interface roughness and sediment volume fluctuations on the resolution of broadband 
pulse propagation in a simple, shallow water environment Predictions of the acoustic 
propagation were accomplished with a broadband, wide-angle parabolic equation model. 
Fluctuations of the bottom interface and volume were computed from realistic spectra 
characterizing each structure. The strength, or level of variability, of each was varied 
independently in order to determine the relative influence on pulse resolution 

T o be consistent with measurement techniques employed in underwater acoustics 
experiments, a beamforming algorithm was written to discriminate vertical arrival angles 
as a function of arrival time. When no fluctuations were introduced, the solution to the 
range-independent problem was found to agree with a standard normal mode analysis 
When fluctuations were introduced, the symmetrical arrival time pattern became distorted, 
indicating a degradation in the resolution of the pulse Such degradations are expected to 
limit the ability of some operational naval systems 

When interface roughness was introduced, the results suggested that the main 
parameter influencing resolution was the ratio of rms roughness to correlation length In 
other words, a high value of rms roughness does not independently imply a degradation in 
the pulse shape if the correlation length is quite large. Similarly, a relatively small rms 
roughness may still significantly affect pulse resolution if accompanied by a very small 
correlation length In contrast, the largest realistic scales of volume fluctuations had little 
impact on the pulse shape A combination of large volume fluctuations and various scales 
of rms interface roughness showed that degradations of pulse resolution are dominated by 
the interface structure. While it may be possible that a different environment (e g , strongly 
downward refracting water column and upward refracting sediment) would increase the 
influence of sediment heterogeneity, such cases would also tend to increase interactions 
with the bottom interface resulting in similar conclusions 
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An attempt was made to quantity these results by measuring the Me - width of the 
pulse at various arrival angles or arrival times as a function of the fluctuation strength. It 
was observed that an rms interface roughness of 10 m with a correlation length of 1000 m 
produced a doubling of the angular width of the arrivals and increased the transmission 
loss by as much as 10 dB The time width was more stable, however In none of the cases 
with volume fluctuations was the signal found to degrade appreciably When the two 
influences were combined, the results were nearly identical to the case with no volume 
fluctuations 

Overall, these results suggest that bottom interactions in shallow water 
propagation are dominated by interface scatter. The degradation of pulse resolution 
depends on a combination of the rms interface roughness and an associated correlation 
length scale. In order to know a priori the limitations of pulse resolution, it is necessary to 
obtain an accurate statistical representation of the sea floor topography. In contrast, the 
sub-bottom volume fluctuations do not appear to significantly influence the forward 
propagation 
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